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Abstract 

Properties of the orbitally excited (L = 1) B® states are studied using l.Ofb -1 
of pp collisions at y/s = 7TeV collected with the LHCb detector. The first ob- 
CN ■ servation of the -B* 2 (5840)° meson decaying to B* + K~ is reported, and the cor- 

responding branching fraction measured relative to the B + K~ decay mode. The 
S s i(5830)° -> B*+K~ decay is observed as well. The width of the £* 2 (5840)° state 
^ ■ is measured for the first time and the masses of the two states are determined 

J-j . with the highest precision to date. The observation of the B* 2 (5840)° ->■ B* + K~ 

decay favours the spin-parity assig nment J p = 2+ for the 5* 2 (5840)° meson. In 
addition, the most precise measurement of the mass difference m{B* + ) — m(B + ) = 
45.01 ± 0.30 (stat) ± 0.23 (syst) MeV/c 2 is obtained. 
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Heavy quark effective theory (HQET) describes mesons with one heavy and one light 
quark where the heavy quark is assumed to have infinite mass pQ . It is an important tool 
for calculating meson properties which may be modified by physics beyond the Standard 
Model, such as CP violation in charm meson decays [2] or the mixing and lifetimes of B 
mesons [3]. It also predicts the properties of excited B and B° s mesons [IH7], and precise 
measurements of these properties are a sensitive test of the validity of the theory. Within 
HQET the B° s mesons are characterised by three quantum numbers: the relative orbital 
angular momentum L of the two quarks, the total angular momentum of the light quark 
j q = |L±2 1 , and the total angular momentum of the B® meson J = |j 9 ±| | . For L = 1 there 
are four different possible (J, j q ) combinations, all with even parity. These are collectively 
termed the orbitally excited states. Such states can decay to B + K~ and/or B* + K~ (the 
inclusion of charge-conjugate states is implied throughout this Letter), depending on their 
quantum numbers and mass values. The two states with j q = 1/2, named £>* and B' sl , 
are expected to decay through an S-wave transition and to have a large (9(100 MeV/c 2 ) 
decay width. In contrast, the two states with j q = 3/2, named -B s i(5830)° and 5* 2 (5840)° 
(henceforth B si and B* 2 for brevity), are expected to decay through a D-wave transition 
and to have a narrow (9(1 MeV/c 2 ) decay width. Table [T] gives an overview of these states. 



Table 1: Summary of the orbitally excited (L = 1) B° s states. 







J p 


Allowed decay mode 
B+K- B* + K~ 


Mass ( MeV/c 2 ) [8] 




1/2 


0+ 


yes 


no 


Unobserved 


B' sl 


1/2 


1+ 


no 


yes 


Unobserved 


B s \ 


3/2 


1+ 


no 


yes 


5829.4 ±0.7 


B* s2 


3/2 


2+ 


yes 


yes 


5839.7 ±0.6 



In this Letter a 1.0 fb _1 sample of data collected by the LHCb detector is used to search 
for the orbitally excited B° s mesons in the mass distribution of B + K~ pairs, where the B + 
mesons are selected in the four decay modes: B + — > J/tp(fi + fi~)K + , B + — > D (K + 7T~)7r + , 
B + — > D°(K + 7r~7T~7T + )7r + , and B + — > D°(K + 7t~)tt + 7t~7t + . Two narrow peaks were 
observed in the B + K~ mass distribution by the CDF collaboration [9]. Putatively, they 
are identified with the states of the j q = 3/2 doublet expected in HQET jl] and are named 
B s i and B* 2 . As the B s i — > B + K~ decay is forbidden, one of the mass peaks observed 
is interpreted as the B s i — > B* + K~ decay followed by B* + — > -B + 7, where the photon is 
not observed. This peak is shifted by the B* + — B + mass difference due to the missing 
momentum of the photon in the B* + — > B +/ ~f decay. While the B* 2 — > B + K~ decay has 
been observed by the DO collaboration as well [10], a confirmation of the B s \ meson is 
still missing. In addition, the B s i and B* 2 quantum numbers have not yet been directly 
determined and their identification in the B + K~ mass spectrum is based on the expected 
mass splitting between the j q = 3/2 states. The B* 2 — > B* + K~ decay has not yet been 
observed but could manifest itself in the B + K~ mass spectrum in a similar fashion to 
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the corresponding B s i meson decay. The B* 2 — > B* + K~ branching fraction relative to 
B* 2 — > B + K~ is predicted to be between 2% and 10%, depending on the B* 2 mass [HhT4"] . 

Recently the Belle collaboration has reported observation of charged bottomonium-like 

states, that could be interpreted as B^B* molecules [15]. To test this interpretation, 
improved measurements of the B* + mass are necessary, and can be obtained from the 
difference in peak positions between B* 2 — > B* + K~ and B* 2 — > B + K~ decays. 

The LHCb detector [16] is a single-arm forward spectrometer covering the pseudo- 
rapidity range 2 < r] < 5, designed for studying particles containing b or c quarks. The 
detector includes a high-precision tracking system consisting of a silicon-strip vertex de- 
tector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift tubes placed downstream. The combined tracking 
system has a momentum resolution (Ap/p), that varies from 0.4% at 5 GeV/c to 0.6% at 
lOOGeV/c, and a decay time resolution of 50 fs. The resolution of the impact parameter, 
the transverse distance of closest approach between the track and a primary interaction, 
is about 20 |xm for tracks with large transverse momentum. The transverse component is 
measured in the plane normal to the beam axis. Charged hadrons are identified using two 
ring-imaging Cherenkov detectors. Photon, electron and hadron candidates are identified 
by a calorimeter system consisting of scintillating-pad and pre-shower detectors, an elec- 
tromagnetic calorimeter, and a hadronic calorimeter. Muons are identified by a system 
composed of alternating layers of iron and multiwire proportional chambers. 

The trigger system [T7] consists of a hardware stage, based on information from the 
calorimeter and muon systems, followed by a software stage that applies a full event 
reconstruction. Events likely to contain a B meson are selected by searching for a dimuon 
vertex detached from the primary interaction or two-, three-, and four-track vertices 
detached from the primary interaction which have high total transverse momentum. These 
are respectively referred to as dimuon and topological triggers. 

The samples of simulated events used in this analysis are based on the Pythia 6.4 
generator [18], with a choice of parameters specifically configured for LHCb [19]. The 
EvtGen package [20] describes the decay of the B mesons, and the Geant4 toolkit 
[2"T||2"2"] is used to simulate the detector response. QED radiative corrections are generated 
with the Photos package [23] . 

In the offline analysis the B mesons are reconstructed using a set of loose selection 
criteria to suppress the majority of the combinatorial backgrounds. The B + — > J/ipK + 
selection requires a B + candidate with a transverse momentum of at least 2 GeV/ c and a 
decay time of at least 0.3 ps. For the other decay modes, the selection explicitly requires 
that the topological trigger, which selected the event, is based exclusively on tracks from 
which the B meson candidate is formed. Additional loose selection requirements are placed 
on variables related to the B meson production and decay such as transverse momentum 
and quality of the track fits for the decay products, detachment of the B + candidate from 
the primary interaction, whether the momentum of the B + candidate points back to the 
primary interaction, and the impact parameter y 2 . The impact parameter y 2 is defined 
as the difference between the y 2 °f the primary vertex reconstructed with and without 



2 
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Figure 1: Invariant mass spectra of the final B + candidates. The signal lineshape is fitted 
with a double Gaussian distribution, while the background is modelled with a second order 
polynomial . (a) _B+ -> J/ipK + , (b) B + -> J D°(A' + 7r-)vr + , (c) 5+ -> _D°(Ji + 7r-7r-7r + )vr + , 
and (d) 5 + — >■ D°(K + 7r~)7r + 7r~7r + decays. The J/^ and _D° masses are constrained to 
their world average values. 

the considered track. 

Following these selections, B + signals are visible above backgrounds in all four decay 
modes. In order to improve their purity, four boosted decision tree classifiers [21] are 
trained on variables common to all four decay modes: the transverse momenta and impact 
parameters of the final state tracks, the transverse momentum and impact parameter of 
the B + candidate, the detachment of the B + candidate from the primary interaction, the 
cosine of the angle between the B + candidate momentum and the direction of flight from 
the primary vertex to the decay vertex, the fit \ 2 °f the tracks, and particle identification 
information. The classifier is trained on data using the sW eights technique [25], with 
the B + candidate mass as a discriminating variable, to unfold the signal and background 
distributions. The cut on the classifier response is chosen by optimizing the significance 
of each B + signal. The final mass distributions for the B + candidates are shown in Fig. [TJ 

The B + candidate mass spectra are fitted using a double Gaussian function for the 
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signal and a second order polynomial for the background. The average mass resolution, 
<7b+, is defined as the weighted average of the Gaussian widths. The B + candidates, within 
a ±2erg+ mass region, are selected for each decay mode. A sample of about 1 000 000 B + 
candidates is obtained and combined with any track of opposite charge that is identified 
as a kaon. 

Multiple pp interactions can occur in LHC bunch crossings. In order to reduce com- 
binatorial backgrounds, the B + and kaon candidates are required to be consistent with 
coming from the same interaction point. The signal purity is improved by a boosted 
decision tree classifier, whose inputs are the B + and the kaon transverse momenta, the 
log-likelihood difference between the kaon and pion hypotheses, and the vertex fit and 
impact parameter \ 2 ■ The training is performed using simulated events for the signal and 
the like-charge B + K + candidates in the data for the background. The same selection is 
subsequently applied to all B + decay modes. The cut on the classifier response is chosen 
by optimizing the significance of the B* 2 — > B + K~ signal. It retains 57% of the signal 
events and rejects 92% of the background events. In order to improve the mass resolution, 
the B + K~ mass fits are performed constraining the J/ip (or D°) and B + particles to their 
respective world average masses [5] and constraining the B + and K~ momenta to point 
to the associated primary vertex. 

Figure |2] shows the mass difference for the selected candidates, summed over all B + 
decay modes. The mass difference is defined as Q = m(B + K~) —m(B + ) —m(K~) where 
m(B + ) and m(K~) are the known masses of the B + and K~ mesons [8], respectively. 
The two narrow peaks at 10 and 67MeV/c 2 are identified as the B s i — y B* + K~ and 
B* 2 — y B + K~ signals, respectively, as previously observed. In addition, a smaller struc- 
ture is seen around 20MeV/c 2 , identified as the previously unobserved B* s2 — y B* + K~ 
decay mode. 

An unbinned fit of the mass difference distribution is performed to extract the Q values 
and event yields of the three peaks. The B* 2 — y B + K~ signal is parameterized by a rela- 
tivistic Breit-Wigner function with natural width V convolved with a Gaussian function 
that accounts for the detector resolution. Its width is fixed to the value obtained from 
simulated events (~ lMeV/c 2 ), increased by 20% to account for differences between the 
B + resolutions in data and simulated events. The lineshapes of the B s i/B* 2 — y B* + K~ 
signals, expected to be Breit-Wigner functions in the B* + K~ mass spectrum, are af- 
fected by the phase space and the angular distribution of the decays as the photon is 
not reconstructed. The resulting shapes can not be properly simulated due to the lack 
of knowledge of the B s i/B* 2 properties. Therefore a Gaussian function is used for each 
B s i/B* 2 —y B* + K~ signals as effective parameterization. The background is modelled by 
a threshold function, f(Q) = Q a e l3 ® +S , where a, (3 and 5 are free parameters in the fit. Its 
analytical form is verified by fitting the like charge B + K + combinations where no signal 
is expected. 

The parameters allowed to vary in the fit are: the yield N b * 2 ^b+k-, the yield 
ratios N Bs1 ^b*+k- /N b ^b+k- and N b ^b*+k- /N b * 2 -+b+k-, the Q values of the 
B sl — y B* + K~ and B* 2 —y B + K~ signals, the mass difference between the B* 2 — y B + K~ 
and B* 2 —y B* + K~ peaks, the natural width of the B* 2 state, the Gaussian widths of 
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m(B + K") - m(B + ) - m(K") [MeV/c 2 ] 

Figure 2: Mass difference distribution m(B + K~) — m(B + ) — m(K~). The three peaks 
are identified as (left) B 3l ->■ B* + K~, (middle) B* s2 ->■ B* + K~, and (right) S* 2 ->■ B + K~. 
The total fit function is shown as a solid blue line, while the shaded red region is the 
spectrum of like-charge B + K + combinations. The inset shows an expanded view of the 
B s i/B* 2 — > B* + K~ signals. The bottom plot shows the fit pulls. 



B s i/B* 2 — > B* + K signals and the parameters of the threshold function. From the yield 
ratios, the relative branching fraction 



B(B* s2 -> B*+K~) N B * 2 ^ B * +K - 



B(B* s2 -)> B+K-) N B: . B+K - 



., ,rel r>B* 



is measured. The B s \ to B* 2 ratio of production cross-sections times the ratio of branching 
fractions of B s \ — > B* + K~ relative to that of B* 2 — > B + K~ is also determined from 

g(pp^B 8l X)B(B 8l ^B*+K-) = N Bsl ^ B , +K - , = „ i/flJj Bsi/B « 2 

a(pp 5: 2 X)i3(^ 2 S+iT-) " ■ 7V S ^ B+ ^ ^ 2 1 J 

These ratios are corrected by the relative selection efficiencies, e r 2 \ = 1.05 ± 0.02 and 
e i 6 2 = 1-03 ± 0.01, using simulated decays. The fit results are given in Table HJ A binned 
X 2 test gives a confidence level of 43% for the fit. 

To determine the significance of the B* 2 — > B* + K~ signal, a similar maximum likeli- 
hood fit is performed, where all parameters of the signal are fixed according to expectation, 
except its yield. The likelihood of this fit is compared to the result of a fit where the yield 
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Table 2: Results of the fit to the mass difference distributions m(B + K )—m(B + )—m(K ). 
The first uncertainties are statistical and the second are systematic. 



Parameter 



Fit result 




fO.46 ± 0.04 ± 0.04 MeV/c 2 
67.06 ± 0.05 ±0.if MeV/c 2 
45.01 ± 0.30 ±0.23 MeV/c 2 
1.56 ± 0.13 ± 0.47 MeV/c 2 

(9.3 ±1.3 ±1.2)% 



(23.2 ±1.4 ±1.3)% 



750 ± 36 
307 ± 46 
3140 ± 100 



of the signal is fixed to zero. The statistical significance of the B* 2 — > B* + K signal is 



A number of systematic uncertainties are considered. For the signal model, the signal 
shape is changed to a double Gaussian function and an alternative threshold function is 
used for the background. The changes in the fit results are assigned as the associated 
uncertainties. The B + decay modes are fitted independently to test for effects that may 
be related to differences in their selection requirements. For each observable quoted in 
Table [21 the difference between the weighted average of these independent fits and the 
global fit is taken as a systematic uncertainty. Additional systematic uncertainties are 
assigned based on the change in the results when varying the selection criteria and the B + 
signal region. The detector resolution of B* 2 — > B + K~ signal is conservatively varied by 
±20%. In addition, the momentum scale in the processing of the data used in this analysis 
is varied within the estimated uncertainty of 0.15%. The corresponding uncertainty on 
the measured masses is assigned as a systematic uncertainty. The uncertainty on the 
determination of the selection efficiency ratios caused by finite samples of simulated events 
is taken as a systematic uncertainty for the branching fractions. Finally simulated events 
are used to estimate the mass shifts of the B s \/B* s2 — > B* + K~ signals from the nominal 
values when the radiated photon is excluded from their reconstructed decays. The absolute 
systematic uncertainties are given in Table [3] while the final results are shown in Table [2J 
The measured B ^b*^b+k-) branching fraction ratio and B* 2 width are in good agreement 
with theoretical predictions ^21414] . 

The mass differences given in Table [2] are translated into absolute masses by adding 
the masses of the B + and kaon [8] and, in the case of the B s i meson, the B* + — B + mass 
difference measured in this Letter. The results are 
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Table 3: Absolute systematic uncertainties for each measurement, which are assumed to 
be independent and are added in quadrature. 



Source 


Q(B sl ) 


Q(B* s2 ) 


m(B*+) - m{B+ 


) t(b; 2 ) 








(MeV/c 2 ) 


( MeV/c 2 ) 


( MeV/c 2 ) 


( MeV/c 2 ) 


(%) 


(%) 


Fit model 


0.00 


0.02 


0.03 


0.01 


0.2 


0.5 


B + decay mode 


0.01 


0.01 


0.02 


0.01 


0.1 


0.1 


Selection 


0.03 


0.02 


0.19 


0.05 


1.1 


0.6 


B + signal region 


0.01 


0.03 


0.11 


0.07 


0.2 


0.4 


Mass resolution 


0.00 


0.01 


0.02 


0.46 


0.2 


0.9 


Momentum scale 


0.02 


0.10 


0.03 








Efficiency ratios 










0.2 


0.2 


Missing photon 


0.01 




0.01 








Total 


0.04 


0.11 


0.23 


0.47 


1.2 


1.3 



m(B* + ) = 5324.26 ± 0.30 ± 0.23 ± 0.17 MeV/c 2 , 
m(B sl ) = 5828.40 ± 0.04 ± 0.04 ± 0.41 MeV/c 2 , 
m(B*J) = 5839.99 ± 0.05 ± 0.11 ± 0.17 MeV/c 2 , 

where the first uncertainty is statistical and the second is systematic. The third uncer- 
tainty corresponds to the uncertainty on the B + mass [8] and, in the case of the B s \ mass 
measurement, the uncertainty on the B* + — B + mass difference measured in this analysis. 

In conclusion, using 1.0 fb" 1 of data collected with the LHCb detector at y/s = 7TeV, 
the decay mode B* 2 — > B* + K~ is observed for the first time and its branching fraction 
measured relative to that of B* 2 — > B + K~ . The observation of the B* 2 meson decaying 
to two pseudoscalars (£>* 2 — > B + K~) and to a vector and a pseudoscalar (B* 2 B* + K~) 
favours the assignment of J p = 2 + for this state. The B* 2 width is measured for the 
first time, while the masses of the B s \ and B* 2 states are measured with the highest 
precision to date and are consistent with previous measurements [3J[TU]- Finally, the 
observed B* 2 — > B* + K~ decay is used to make the most precise measurement to date of 
the B* + — B + mass difference. This measurement, unlike others reported in the literature, 
does not require the reconstruction of the soft photon from B* + decays and therefore 
has significantly smaller systematic uncertainty. High precision measurements of the B* + 
mass are important for the understanding of the exotic states recently observed |15j . 
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